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Abstract — A modular fluid-flow model for network congestion 
analysis and control is proposed. The model is derived from 
an information conservation law stating that the information 
is either in transit, lost or received. Mathematical models 
of network elements such as queues, users, and transmission 
channels, and network description variables, including send- 
ing/acknowledgement rates and delays, are inferred from this law 
and obtained by applying this principle locally. The modularity 
of the devised model makes it sufficiently generic to describe 
any network topology, and appealing for building simulators. 
Previous models in the literature are often not capable of 
capturing the transient behavior of the network precisely, making 
the resulting analysis inaccurate in practice. Those models can 
be recovered from exact reduction or approximation of this 
new model. An important aspect of this particular modeling 
approach is the introduction of new tight building blocks that 
implement mechanisms ignored by the existing ones, notably 
at the queue and user levels. Comparisons with packet-level 
simulations corroborate the proposed model. 

Index Terms — Congestion control modeling;Fluid-flow models; 
Queueing model; Self-clocking; 

I. Introduction 

Network modeling is challenging due to the very hetero- 
geneous nature of communication networks, mixing physics, 
electronics and computer science. This heterogeneity coupled 
with intrinsic properties of physical and mathematical laws 
prevent the development of an efficient bottom-up approach. 
This is why finding macroscopic laws capturing critical phe- 
nomena is of interest. These laws should provide an abstrac- 
tion of the microscopic level by identifying and relating the 
fundamental macroscopic network parameters. 

We derive a modular network model using three fundamen- 
tal laws. The first law is a packet conservation law that facili- 
tates the derivation of models as building blocks and simplifies 
their mathematical expression. The second law defines a model 
for queues and, finally, the last one concerns the existence of 
a user model. We will show that each law has implications 
for the network modeling problem and, more importantly, 
will allow to solve yet unresolved problems, especially at 
user level. One important property of the developed model 
is its modularity. Indeed, the modeling technique allows to 
develop each element independently of the others, leading then 
to building blocks which may be interconnected as desired 
to build any network topology. The model shares exactly 
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the same structure with a real network which is built by 
interconnection of several elements. This property is also very 
appealing for simulation purposes where a topology can be 
easily simulated by connecting the building-blocks. 

The proposed model is developed in several steps. The first 
one is the modeling of lossless transmission channels directly 
from the first law. The second step is the derivation, again 
from the first law, of the so-called ACK-clocking model (TJ, 
which ties flow, flight-size and round-trip time (RTT) together. 
This result is of great importance in network modeling. 

Based on the first and second laws, a causal RTT expression 
is developed in the third step. This causal RTT expression 
however requires an extension of the buffer model. Indeed, 
there are two main limitations to the buffer model usually 
considered in the literature (and as stated in the second law). 
First, the model does not explicitly define the queue as a FIFO 
queue (i.e. order preserving) in which the packets maintain 
their relative positions. An internal buffer description should 
capture this, at a flow level. To this aim, the flows should 
therefore be considered as very viscous repelling liquids which 
do not mix. The second limitation concerns the solving of the 
output flow separation problem, primordial for the description 
of buffer interconnections and, as we shall see later, for the 
derivation of an exact expression for the acknowledgment 
flows. 

Finally, the last step is devoted to the derivation of a 
complete user model, based on the first and third laws. This 
part constitutes one important contributions of the paper. 
Indeed, the conversion of congestion window size into flow has 
been a major obstacle preventing the improvement of network 
models. The static-link model Q assimilates the flows to be 
equal to the derivative of congestion window sizes. It has good 
modeling properties in the absence of cross-traffic and when 
propagation delays are homogeneous. It is however rather 
inaccurate in more realistic scenarios. This validity domain is 
theoretically proved in this paper by showing that the proposed 
model reduces to the static-link one when some conditions are 
met. The integrator link model Q, JU, 0, J6], Q improves 
the description by correcting the irrelevant behavior of the 
model when affected by cross-traffic. Yet, some characteristics 
of the buffer response were not well captured: the response 
speed and the high slope when the congestion window size 
increases. The joint-link model O, J9) consisting of merging 
the static-link and integrator-link models has been introduced. 
This approach improves the network model by capturing some 
characteristics unmodeled by the previous descriptions. It has 
been shown that these flow models can, in fact, be considered 



as approximations of the ACK-clocking model (8), (9), ifTUl 
from which higher order approximations can also be defined. 
More recently, the ACK-clocking model has been exactljfl 
considered in [H, IfTUl and has lead to important improvements 
in terms of precision. In this paper, we also do not make 
any approximations and use the ACK-clocking model in a 
new fashion, leading to a new user model explicitly using the 
received flow of acknowledgements. This new model exactly 
captures both the ACK-clocking and the decreasing of the 
congestion window size, the latter being not captured by the 
existing models. 

The modeling technique is applied to a single- 
buffer/multiple-users topology for which it is possible 
to show that, under some certain conditions, the static-link 
model can be naturally recovered. This shows that the 
static-link model is more general than it was previously 
known Q, iflOl . Unlike other versions of the static-link 
model, the obtained one involves a state-dependent time-delay 
ifTTl representing the queuing delay lfl2l . 

Finally, simulation results are provided for some simple 
topologies to compare the proposed model with previous 
ones. It is shown that the model reproduces the results of 
0, IflOl . as expected since those models are based on the 
ACK-clocking model. By transitivity, the proposed model also 
matches NS-2 simulations and the experiments reported in 
|8|, 1 10]. Simulations also illustrate that the proposed model 
is also able to capture the behavior of the queues and users 
when the congestion window size decreases, as opposed to the 
existing ones. Another important property of the model is its 
modularity which makes it easy to manipulate when building 
network topologies and setting up simulations, where blocks 
just have to be connected to others. 
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II. Definitions and laws 

A. Networks and Graphs 

It is convenient to introduce here the particular network 
graph representation considered in the paper. It is different 
from the regular ones since it places all network elements on 
graph edges, leaving nodes with the role of connecting points, 
as in electrical circuits. Four types of nodes are distinguished: 



+ 



bj for user 



the input nodes , bj and output nodes 
Ui and buffer bj, respectively. The superscripts have to be 
understood as a temporal order of reaction or causality: the 
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Fig. 3. Example of graph with 4 edges: one user edge u = (u~,u + ), 
one buffer edge b = and two transmission edges l\ = 
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data come at (-) and leave at (+). We denote any edge E of 
the graph by (x, y) where x and y are the input and output 
nodes respectively. Moreover, given any edge E, the input and 
output nodes are given by (3(E) and e(E) respectively. 
According to these definitions, a queue edge is always 

, uf) and a trans- 



denoted by (b i , b~l 
mission edge by (b^ 



a user edge by 



if,bj) or (bf,b^),i^k. This 



is illustrated in Fig. [3] We call a circuit, say C, a path from 
the output to the input of a user, i.e. C = (u + , u~). In Fig. [3] 
the only possible circuit is given by C = (u + , b~ ,b + ,u~). 

B. Fluid-Flow Paradigm 

In this paper, we shall present several laws allowing for 
the congestion modeling in a fluid-flow framework lfl3ll . 0, 
E), 0. Using such a groundwork, the modeling problem 
admits a scalable solution in which the network asynchrony 
is captured through appropriate expansions and compressions 
of the time-space. It also enables the use of well-established 
mathematical tools: real functions analysis, integration theory, 
dynamical systems, difference and differential equations, etc. 
A first conclusion is that a proper definition for flows of data 
can easily be stated. 

Definition 2.1: Given any edge E of the network, a non- 
negative integrable scalar flow (j>(x, s), (x, s) <G E x R + and 
time instants to, t 6 R+, t > to, we have 



N x (t,t ) 



4>{x, s)ds 



(1) 



where N x (t, to) is a packet counter, i.e. the number of packets 
that have passed through point x between to and t is given by 
N x {t, to), and the integral is a standard one, e.g. the Lebesgue 
integral. 

Unlike electrical circuits where one flow (the current) cir- 
culates on the edges, a vector-valued flow circulates on the 
edges of the graph representation of the network. Using the 
notation defined in Section III- Al we can build the flow vectors 



(x,t) using the 'col' operato 



1; 



(x,t) = coi^ y> k (x,t)] 



(2) 



Although the definitions for the ACK-clocking model slightly differ. 



2 which stacks in column the arguments. 



where u(x) is the number of flows in parallel using node 
x. Note that the nodes then implement multiplexers and 
demultiplexers for appropriate routing of the flows. 

C. Law 1: Information Conservation Law 

The first law is essentially a conservation law relating flow 
integration on two different domains. This conservation law 
follows from the remark that the quantity of information is 
preserved in a communication network: the data is either in 
transit, lost or received. Assuming lossless networks, it is 
possible to determine the total number of packets in transit on 
any edge, simply by counting the number of entering packets 
over a certain time horizon. 

Law 1: Given any edge E of the network, then for all t G 
R+ there exists a time to(t) G K+, to(t) < t such that 

P E (t) := L<t>(o,t)de 

= f to{t) ct>{P{E),s)ds (3) 

= N 0(E) (t,t o (t)) 
The integration over E is an abstract integral which has to be 
understood as a flow integration from /3(E) to e(E), that is, 
the number of packets P E (f) in the edge E = (/3(E), e(E)) 
at time t. 

The main features of this law are the domain of integration 
exchange and the discretization of the spatial domain to 
nodes only. These considerations dramatically simplifies the 
modeling since it is no longer necessary to consider the flows 
at any point x G E but only at input nodes /3(E). This is 
illustrated by the following proposition: 

Proposition 2.1: The input flow cf>(/3(E), ■) and the output 
flow cf)(e(E), ■) of edge E verify 



4>(e(E),t)=t (t)'t(f3(E),to(t)). 



(4) 



Proof: Since Np( E ^(t,to(t)) is me current number of 
packets in the edge E at time t, then differentiation with 
respect to time provides the corresponding rate of variation 

[N m) (*, t (t))}' = 4>(P(E),t) - t (t)'<i>(p(E), t (t)). 

Moreover, the variation of the number of packets verifies 

[N p{E) (t,t (t))}' = cf>(/3(E),t) - <]>(e(E),t) 

which is nothing else but the difference between the input 
and output flows. The result follows from identification of the 
equalities. ■ 
This proposition turns out to be very useful to derive models 
for transmission channels and buffers. 

D. Law 2: Queues are Flow Integrators 

The law given below defines the behavior of queues in- 
volved, for instance, inside routers and servers. Following past 
works and our understanding of the problem, the integrator 
model for queues is the most realistic. 

Law 2: The queue dynamics of buffer i is governed by the 
model 



with aggregated output flow rate 



if Ci(t) 
otherwise. 



(6) 



Above, qi, Ci and 4>j(b~,t) represent the queue size, the 
maximal output capacity and the flow of type j at the input, 
respectively. The condition Ci(t) is given by 



d(t) := [q t (t) > 0] V 



-,t) > a 



(7) 



Qi(t) 



<f>j{K,t)-n(t) 



(5) 



The corresponding queuing delay can be easily deduced using 
the relation Ti(t) = qi(t)/ci. 

The above model can be proved to be the small packet limit 
of an M/M/l queue (9). It can also be refined to capture 
additional features like finite maximal queue length, flow 
priorities, multiple output capacities, etc. 

E. Law 3: Users Model Existence 

The last law concerns the user protocol description and the 
way it dynamically reacts to congestion in the network. 

Law 3: There exist bounded functions Vi, Wi and Ui such 
that the trajectories (zi(t),Wi(t)) of the following continuous- 
time model defined over t £ M+ 

Wi (t) = m(zi(t),m(t)) (8) 

cf>i(uf,t) = Ui(wi(t),(j}i(uJ ,t)) 

match the trajectories of the asynchronous protocol at points 
in R + fl Tj. Above, Z{, (f>i(u^,-) and 4>i(u^ 7 -) are the 
state of the protocol, the measurements, the acknowledgment 
flow rate and the user sending flow respectively. The con- 
gestion window size Wi is considered here as the number of 
outstanding packets to track and is supposed to be (weakly) 
differentiable. 

A procedure to solve the above interpolation problem has 
been first proposed in (4) for TCP and reused in |14] Appendix 
C] for FAST-TCP. 

III. Preliminary results 
A. Transmission Channel With Constant Propagation Delay 

The following result concerning transmission channels is an 
immediate consequence of law Q] 

Result 3.1: Given a lossless transmission channel, corre- 
sponding to an edge E, with constant propagation delay 
T > 0, the output flow is given by 

<f>(e{E),t) = <i>{P(E),t-T). (9) 

Proof: Following law Q] the number of packets in transit 
P E (t) in the edge E at time t G M + obeys 

P E (t) = L^(x,t)dx 

= r tQ{t) 4>(p(E),s)ds (10) 
- Np (E) (t,t (t)), 

3 The asynchronous discrete decision instants of the actual protocol are 
assumed to belong to a countable set Tj. 
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where to (i) = t — T since the propagation delay is constant. 
Indeed, a packet sent at time t — T will, at time t, still be 
on the edge but about to leave. The result follows then from 
Proposition 12.11 ■ 

B. The ACK-Clocking model 

The ACK-clocking model 0, fT|, HO) is certainly the most 
important consequence of law Q] This model characterizes the 
flight-siz^j F i(Ci,t) := Pd{t) of a user i at any time t £ 
R + over a closed circuit C, L = (u^,u~). The importance of 
the ACK-clocking model lies in the semantic it adds to the 
model by relating RTT, flow and flight-size togethetQ. From 
this result, it will be possible to derive a number of important 
properties and rules for the users. 

Result 3.2 (ACK-Clocking): The ACK clocking model is 
given by 



Fi{Ci,t + RTTj{£}) = / <t>(0, t + BHTi{t})d9 

JCi 
t+RTTi{t} 



<j>i(uf,s)ds 



(11) 

where RTTijt} is the RTT of a packet sent at time t in the 
circuit Ci by user i, A 
Proof: We assume here that the circuit is lossless for 
simplicity. The flight-size is indeed a spatial integration of 
flows since the number of packets in transit is equal to the 
spatial integral of the flows over the circuit. Thus, following 
law [1] it is possible to convert the spatial integration into a 
temporal one provided that we can determine the integration 
domain. To obtain it, we use the notion of RTT and suppose 
that a data sent by user i in the circuit Cj at time t has a 
round-trip-time given by RTTj{t}. This means that the data 
sent between t and t + RTT i{t} are still unacknowledged at 
t+RTTi{t}, thus still in the circuit. Hence, the corresponding 
temporal integration has bounds t and t + RTTi{t}. ■ 

IV. RTT EXPRESSION AND INTERNAL FIFO BUFFER 
MODEL 

In the light of the discussion above, it turns out that a 
model for the RTT is necessary in order to characterize and 
compute the flight size F%(Ci,t + RTTi{t}). A first step 

4 The number of outstanding packets. 

5 Note however than in 1101 the congestion window size is considered 
instead of the flight-size, which is rather different. Equivalence holds when 
some conditions, such as the condition flight-size>congestion window size, 
are met. We will come back on this in Section IV Bl 



forward towards a RTT expression is the analysis of queuing 
delays in buffers. 

A. Forward/Backward Delays and Causal RTT Expression 

Since RTT directly depends on the queuing delays, its 
computation essentially relies on calculating the queueing 
delays. The results are taken from a previous work of us lfT2l 
with the difference that we relate them here to law [2] 

Considering the buffer model (O with queueing delay Ti(t), 
we define the forward delay operator fi(t) := t + t,(£) 
which maps, at a flow level, any packet input time t to the 
corresponding packet output time fi(t). Although it is easy 
to derive and understand, this operator leads to a noncausal 
RTT expression which is not desirable. To observe this, let us 
consider a closed circuit C with N queues, indexed from 1 to 
TV. The indices and N + 1 are used to denote the input and 
output of the circuit respectively. Given a packet input time t, 
the corresponding packet output time tc is given by: 



tc(t) = &c(t) 



= R-nn+i o f N oTZ 



N-1,N 



/iV-1 ° 



o7 ^2,3 ° H ° ^1.2 ° h ° ^0,1 



(12) 



where TZ^j corresponds to the transmission channel operator 
(constant delay operator) between queues i and j, ev is the 
evaluation map and o is the composition operator. The same 
formula, albeit expressed in different ways, has been also 
obtained in |H1 Section 3.3.5] and EoJ equations (7d-7f)]. 

It is clear that &c is noncausal since it requires the 
knowledge of future information, which is not available. In 
order to solve this problem, the backward delay operator 
|[T2l expressing the input time t as a function of the output 
time tc, is considered instead. Hence, the problem reduces 
to inverting the forward delay operator. The existence of this 
inverse operator and some of its properties are recalled below: 

Result 4.1 <H2^): The operator fa is invertible if and only 
if the input flow of the corresponding buffer is almost every- 
where positive. A 

Result 4.2 f/[72V): The functions gi := f7~ obeys: 

9i(t) = t-Ti(gi(t)) 

g'iit) 01 I X] tkiK 1 9i(t)) ] i£Ci(gi(t)) 



fc=i 



1 



otherwise 



where g[(t) stands for the the upper right Dini derivative of 
9i {t), i.e. D+[g,](t) = limsup H0 h' 1 (g t (t + h) - g t (t)) . A 
Using the backward delay operators gi, the packet sending 
time t can be computed from the reception time tc through 
the causal expression: 

t(t c ) = ®c{tc) 

SS C = ^0,1 ° 9x ° ° 92 o 72-2,3 ° •• • (14) 
°9n-i ° 7Ln-i,n ° 9n o Hn,n+i ° ev. 



Example 4.1: In the single-user/single-buffer case, the ex- 
pressions reduce to 

t(t c ) = g(t c -T b )-T f 

t c = t + Tf+T b + T(t + T f ) 

RTT{<} = t c -t (15) 

= tc-g(t c -T b )+T f 

= T b +T f +T(g(t c -T b )) 

where Tf and T b are the forward and backward propagation 
delays corresponding to 7?-o,i and IZ-i^- 

Using the backward expression of the RTT, it easy to obtain 
the following result: 

Result 4.3: The flight size obeys 

(t) 

Fi{Ci,&cM) = J(f>(u+, S )ds 

(16) 

Fi(C l7 t) = J <f>{u+,s)ds. V 

m Ci {t) 

This is a direct consequence of law Q] (through the ACK- 
clocking model) and law [2] However, the problem is not 
completely resolved yet since it is still rather unclear how the 
queuing delays along a given circuit may be computed. Indeed, 
calculating the queueing delays requires the knowledge of 
the input flows of all buffers, and hence necessitates a way 
of splitting the upstream buffer aggregated output flows into 
distinct 'atomic' flows. Otherwise, a modular description of 
buffer interconnections is not possible. 

B. FIFO Buffer Output Flow Separation 

Without further considerations on the queue type, there 
exists an infinite number of ways to separate the aggregated 
output flow directly from the queuing model of law [2] When 
a FIFO queue (i.e. order preserving) is considered, it turns 
out that the output flow separation problem is easily solvable 
using laws Q] and [2] The FIFO characterization and output flow 
separation problems have been fully solved in ifPH . In this 
section, we will recall and explain these results and connect 
them to laws Q] and [2] 

Result 4.4 (472V): Let us consider the queueing model (fSJ 
which we assume to represent a FIFO queue. The output flow 
corresponding to the input flow 4>g(b~ ,t) is given by 



Forward queueing delay r 



MK,t) 



if Ci(t) 



otherwise 



(17) 

Proof: The proof is available in |[L2l and is based on the 
analysis of the contribution of each input flow to the queue 
size. Then, using laws [T] and |2] it is possible to split the 
aggregate output flow into atomic output flows that correspond 
to each input flow. ■ 
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Fig. 5. Queue/Buffer block 



This modef] deserves interpretation: the output flows consist 
of scaling and shifting of the input flows. The delay accounts 
for the high flow viscosity and captures the queue FIFO 
behavior, at the flow level. This model also tells that the output 
flow corresponding to the input flow 4>i(b~ ,t) is expressed as 
a (delayed) ratio of the input flow <fie(b~,t) to the total input 
flow that entered the buffer at the same time. Hence, the output 
flows are proportional to relative flows modeling the probabil- 
ity of having a packet of certain type served at time t. This 
probability is then scaled-up by the maximal output capacity 
to utilize the available bandwidth. This nonlinear expression 
for the output flows describes the flow-coupling, through the 
nonlinear expression in dl7l . and clock-coupling phenomena 
discussed in Section III-BI through the delay dynamical model 
depending on input flows. Note also that the complete model 
(|5]>-(fT7]i is a quite unusual dynamical system consisting of a 
delayed direct feedthrough operator whose delay is implicitly 
defined by a dynamical expression. 

Example 4.2: To illustrate the model, let us consider two 
input flows given by <fri(t) = 0.55c (1 + Sq(£)) and <p2(t) = 
0.55c (1 - Sq(i))where c = lOOMb/s and the function Sq(i) 
is a square function of period T = Is. Since the flows are 
in phase opposition and their sum constant, they lead to an 
alternation of packet types in the queue while respective packet 
populations remain very close to each other at any time. 
Therefore, the output flows should reflect the contents of the 
queues and the model should be able to keep track of the order 
of arrival of packets in the queue. The simulation results are 
plotted in Fig. [6] where we can see that the proposed model 
reflects exactly the contents of the queue. Indeed, the output 
flows are deformations of the input flows: they are scaled 
in amplitude and time, so that the surfaces below the curves 
(the number of packets) over one period are the same. This 
illustrates the law of conservation of the information. In Fig.|7j 
a comparison of the results obtained by the model d5ll-([T7li and 
NS-2 is shown. To overcome the problem of flow computation 
in NS-2, the number of input and output packets of the queue 
is computed for both users and compared. We can see that 
the model matches perfectly NS-2 simulations showing that 
the model well captures the internal organization of packets 
inside the queue. This also means that the output flow model 
is relevant and represents what really happens at the output of 
the queue. 

6 Such a model has also been proposed in 1151 . 1161 without any proof. It is 
proved here that this model is an immediate consequence of the conservation 
law [T] 




Fig. 6. Predicted output flows. Plain: Fig. 7. Model predicted and NS- 
output flows, dashed: input flows 2 simulated input and output number 

of packets: flow 1 (top) and flow 2 

(bottom). 



V. A New user model 

As explained in the introduction, the user modeling problem 
is partially an open question and a solution, based on Laws 
Q] and [3] is proposed here. The distinct notions of flight- 
size, ACK-flow and congestion window size are clarified 
first and associated with each other. Then, the problem of 
computing users sending flow <p(uf,t) is solved. Finally, 
the congestion-window-to-flight-size conversion problem, ac- 
counting for flight-size rate of variation constraints, is ad- 
dressed. 
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Fig. 8. User block 



A. ACK-Clocking Dynamics and User Flow Computation 

The approach exposed here is based on the ACK-clocking 
model of Section IIII-BI obtained from Laws [T] and [2] 

Result 5.1: Let us consider a circuit d = (uf,u^). Then 
the ACK-flow the user m receives is given by 



{ur,t) = 3g' Ci {t)<t> i {<^cM)- 



(18) 



Proof: The key idea is to remark that Fi(Ci,t) = 
N +(t, 3§Ci(t))- Hence, using Proposition 12.11 and noting that 
the ACK-flow is the leaving flow from the circuit, we obtain 
the result. ■ 
Note that the differentiation of ( [Tol l also yields 

4> i {uf,t) = F' i {C i ,t)+^' Ci {t) ( j> i {ut,^cM (19) 

meaning that to maintain the same flight size the user has to 
naturally send data at the same rate as receiving ACK packets: 
this is exactly ACK-clocking but expressed at the flow level. 
By flow integration, we can easily recover the 'packet-level 
ACK-clocking'. A similar formula, albeit not derived from 
any conservation law, is considered in ifTTl . It however uses 
a different model than <JTVj for representing the buffer output 
flowfO and involves the congestion window size instead of the 

7 It can be shown that the output flow model of 1171 does not match NS-2 
simulations in general, but only in some specific cases, such as the constant 
flows case. 



flight-size. 

B. User Flow, Flight-Size and Congestion Window Size 

We clarify here the relation between the user sending flow 
(j)(uf,t) and the congestion window size Wi(t). First, recall 
that the congestion window size corresponds to the desired 
flight-size while the flight-size is the current number of packets 
in transit. The congestion window size is then a reference to 
track while the flight size is the controlled output. The control 
input is the user sending flow for which constraints must be 
considered. 

Indeed, when the congestion window size increases, the 
user can immediately send a burst of new packets to equalize 
the flight-size and congestion window size. In such a case, 
we can ideally assimilate them to be equal (and so are their 
derivatives). The small delay corresponding to the protocol 
reaction time can be easily incorporated in the constant part of 
the RTT The problem is, however, slightly more difficult when 
the congestion window size decreases and becomes smaller 
than the flight-size. In such a case, we cannot withdraw packets 
from the network and the only thing we can do is to wait 
for new ACK packets until the flight size becomes equal to 
the congestion window size. Therefore, while the slope of the 
flight-size is unbounded from above, it is basically bounded 
from below. In (8), a rate-limiter is used to control the slope 
of the flight size but this is rather limited due to the absence 
of any ACK-flow model and the difficulty arising from the 
time-varying nature of the lower bound value to consider. In 
most recent works, such as IfTTl . |[T0l . this problem is excluded 
by considering that the flight-size is always smaller than the 
congestion window size and that the congestion window size 
does not decrease 'too rapidly'. We provide here an explicit 
solution to this problem based on a hybrid modeling of the 
user behavior. 

According to the above discussion, the flight-size must obey 



Fi(Ci,t)' = 



Witt) 

(u 



if Ti(t) 
otherwise 



(20) 



where %{t) is a condition which is false when the flight-size 
must be decreased and true otherwise. 

Result 5.2: The flight-size F i(Ci, t) satisfies ( f20l > if the user 
sending rate is defined as 



Wi(t) + ,t) 





if Ti(t) 
otherwise 



(21) 



where %{t) = ([^(t) = 0] A [w t (t) + <t>{ur ,t) > 0]) and 





w t {t) + 4>{u- ,t) 



if Ti(t) 
otherwise. 



(22) 



Moreover, this model is the simplest one. A 
Proof: The virtual buffer 7Tj, taking nonpositive values, 
measures the number of ACK packets to retain in order to 
balance the flight-size and congestion window size. When the 
virtual buffer has negative state, i.e. iri(t) < 0, the arriving 
ACK-packets have to be absorbed until the state reaches 0. 
Once zero is reached, the user can start sending again until 



the congestion window size decreases too rapidly, i.e. Wi(t) < 
—<fi(u~,t). Substitution of the user sending rate defined by 
d2"fl-ll2"2l in (O yields the flight-size behavior $2&. To see 
that the model is minimal, it is enough to remark that both 
conditions in 71 (t) are necessary. ■ 
The state of the extended user model hence consists of both 
the state of the congestion controller Zi and the ACK-buffer 



VI. Example of a single buffer topology 

A. The Single-Buffer/Multiple-User Topology with FAST-TCP 
protocol 

In this section, we consider a single-buffer/multiple-users 
topology interconnected by lossless transmission channels. 
The forward and backward propagation delays of user i are 
denoted by T/ and T b respectively. For illustration, we use 
the following FAST- TCP model as user protocol model 



Wi(t) = 7 



r(9 l (t)) 



v l (t) 



(23) 



reduces to 

f(t) 

Fi(t) 

0K + ,*) = 
</>K~,t) = 
v(t) 

9b(t) 



c -1 7/(*) - 1 



Wi{t) + <j)(ur,t) 
J t gb{t) _ Tf ^ui,9)d9 

r wi(t) + 4>(u] ,t) 

I o 

( c<P(u+ g b (t)~Tf) 
\ T,j 4>(u+,g b (t)-Tf) 

{ ( j ) (ut,t-T b -Tf) 
9(t T b ). 



if C(t) 
otherwise 
if %{t) 
otherwise. 



if 75(f) 
otherwise 

if C{t) 
otherwise 



(25) 

Assuming the buffer is always congested (i.e. C{t) holds true) 
and all the users are active (i.e. the 7I(t)'s are true) we obtain 



T{t)=c- 1 Y J W l {t-Tf) 



(26) 



which is the static-link model. Integrating the above equation 
from to t we obtain 



r(f) 



(27) 



where Wl (t), T 4 = T( + T? and g^t) = g(t - T?) are the 
congestion window size, the total propagation delay and the 
backward queuing delay respectively. 

1 ) General Model: The general model is given by (|23l and 



f(t) 
Fi(t) 

0K + , 



t) = 




if C{t) 
otherwise 
if Ti(t) 
otherwise. 



if Ti(t) 
otherwise 



if C(t) 
otherwise 



(j>(u t ,t) 

v(t) 

(24) 

where S(t) denotes the normalized cross-traffic 8(t) S [0, 1). 

2) Homogeneous Delays and No-Cross Traffic - The Static- 
Link Model: In ITOl . it is shown that the ratio-link and the joint 
link models are actually approximations of the ACK-clocking 
modeO, i.e. 0th and 1st order approximations of it. The static- 
link model |2| is however retrieved after linearization and Pade 
approximations of exponentials corresponding to delays. This 
'proof makes the static -link model only valid very locally. 
Below we show that the static link model is exact in some 
certain situations and we give a proof for its domain of validity. 

In the case of homogeneous delays, i.e. TV = , T* = T b , 
i = 1, . . . , N, and absence of cross-traffic, i.e. 5 = 0, model 



where we assumed r(0) = and Wi(0) = 0, i = 1, . . . , N. 
Thus, according to the proposed model (l24l . the static-link 
model is valid whenever the necessary and sufficient condi- 
tions are met 

• the buffer is permanently congested, i.e. C(t) holds true; 

• the propagation delays are homogeneous, i.e. t( = T?, 
T b = T b , i = 1,...,N; 

• the cross-traffic is absent, i.e. (5 = 0; 

• the users are not in ACK-retaining mode, i.e. 77 (^) holds 
true for all i = 1, . . . , N. 

By substituting the above static-link model ( 1271 ) into the 
user model (l23l with homogeneous delays, we obtain a new 
model involving the dynamics of the window sizes only. The 
local stability of this model has been studied in detail in IfTlfl . 

VII. Model validation 

The considered scenarios for validating the proposed model 
are taken from (8), iflOl . The results from NS-2 have been 
slightly shifted in time so that the congestion window variation 
times match. Unlike ifTOl . the results are not shifted vertically 
and hence, this results in small discrepancies. If, however, 
the NS-2 results were shifted vertically so that they match 
initial equilibrium values, then the curves would match almost 
perfectly. 

A. Single-Buffer/Multiple-Users 



6+ 





involving the congestion window size instead of the flight-size. 



Fig. 9. Topology of scenarios 1 & 2 



In this section, we consider the interconnection of two users 
through a single resource, as depicted in Fig. [9] The bottleneck 
has capacity c = lOOMb/s and the packet size including 
headers is 1590 bytes. The following scenarios from (8), ifTUl 
are considered: 

• Scenario 1: the congestion window sizes are initially 
iv® = 50 and w§ = 550 packets, at 3s w\ is increased to 
150 packets. The propagation delays are T\ = 3.2ms and 
T2 = 117ms for users 1 and 2 respectively; see Fig. [10] 

• Scenario 2: the congestion window sizes are initially 
iv® = 210 and w° = 750 packets, at 5s w\ is increased 
to 300 packets. The propagation delays are T\ = 10ms 
and T2 = 90ms for users 1 and 2 respectively; see Fig. 





Fig. 13. Scenario 3: queue 1 (top) Fig. 14. Scenario 4: queue 1 (top) 
and queue 2 (bottom) and queue 2 (bottom) 

Again, the model is able to capture the network behavior well 
and retrieve previous results given in J8), IfTUl . This is due 
to the equivalence between the ACK-clocking models in this 
case. 

We introduce now a constant cross-traffic x c \ = Ci/2 on the 
first link. Initiall>0, we set iu? = 1200, wl = 1600, w° 3 = 5 
and we consider the following scenarios: 

> Scenario 5: The congestion window wi is increased by 
200 packets at 10s; see Fig. Q3] 

> Scenario 6: The congestion window W2 is increased by 
200 packets at 10s; see Fig. [161 



Fig. 10. Scenario 1: Queue size Fig. 11. Scenario 2: Queue size 

We can see that the obtained results match quite well NS- 
2 simulations and are very close to the results reported in 
JH, ifTOl . This is expected since both models are based on 
conservation laws which are identical when the flight-size is 
smaller than the congestion window size. Note also that the 
model captures both the transient phase and the steady-state, 
unlike the static-, ratio- and joint-flow models. 

B. Multiple-Buffers/Multiple- Users 




Fig. 12. Topology of scenarios 3 to 6 (8 represents cross-traffic I/O nodes). 

Here, we consider the case of two buffers interconnected 
in series (see Fig. [T2l with capacities c\ = 72Mb/s and C2 = 
180Mb/s. The packet size including headers is 1448 bytes. 
The link propagation delays are 20ms for link 1 and 40ms 
for link 2. The total round-trip propagation delays are T\ = 
120ms, T2 = 80ms and T3 = 40ms for sources 1, 2 and 3 
respectively. Initially, the congestion window sizes are w\ = 
1600 packets, w§ = 1200 packets and = 5 packets. The 
following scenarios from (8), IfTUl are considered: 

• Scenario 3: No cross-traffic and the congestion window 
wi is increased by 200 packets at 10s; see Fig. Qj] 

> Scenario 4: No cross-traffic and the congestion window 
W2 is increased by 200 packets at 10s; see Fig. [141 




Fig. 15. Scenario 5: queue 1 (top) Fig. 16. Scenario 6: queue 1 (top) 
and queue 2 (bottom) and queue 2 (bottom) 

Again the obtained results are identical to the results given 
in 0, IfTUl and match NS-2 simulations. Notice the reaction 
time, consisting of the sum of the queueing delay and the 
propagation delay of link 1, between the moment when the 
congestion window size increases and the moment when the 
second queue sees the flow increases. This illustrates that the 
model captures well the communication path and the order of 
elements. 

C. Decrease in Congestion Window Size 

The models proposed in ifTTl . IfTUl do not capture sudden 
decreases in the congestion window size that would cause 
the buffer to empty or become smaller than the actual flight- 
size, that is, smaller than the number of packets in flight. The 
proposed model captures these phenomena since 1) the ACK- 
clocking model derived from the conservation law involves the 
flight-size rather than the congestion window size iflOl ; and 2) 
the user model involves an ACK-buffer to count the number 
of packets to remove from the network before starting to send 
again. Note that the derivation of the user model including 

This scenario is actually identical to the one in flOl Section III.B.3], the 
congestion window sizes initial values in the paper are inexact. 



the ACK-buffer has been possible thanks to the availability 
of an explicit expression for the flow of acknowledgments, 
itself made computable through an explicit solution for the 
queuing delay and the buffer output flows lfl2l . ifTHl . In (8), 
this case is handled by adding a rate limiter to constrain the 
(negative) slope of the queue size. It is however rather difficult 
to characterize due to the time-varying nature of the lower- 
bound on the slope which depends on the received rate of 
acknowledgment, unfortunately unavailable in the framework 
of the thesis 10. 

Let us consider the single-user/single-buffer case where 
the total propagation delay is T = 150ms, the packet size 
including headers is 1040 bytes and the initial value of the 
congestion window size is w° = 500. A t = 5 seconds, the 
congestion window size is halved. We consider the following 
scenarios 

• Scenario 7: c = 12.5Mb/s and no cross-traffic; see Fig. 

Oil 

• Scenario 8: c = 25Mb/s and half capacity used by cross- 
traffic; see Fig. [181 




Fig. 17. Scenario 7: Queue size (top) Fig. 18. Scenario 8: Queue size (top) 
ACK buffer (bottom) ACK buffer (bottom) 

We can see that we obtain exactly the same results as NS- 
2 simulations (and the rate-limiter model reported in (8)). 
As desired, the ACK-buffer measures (counts) the number of 
packets to remove before starting to send again. 

VIII. Conclusion 

This paper presents a modular fluid-flow model to analyze 
congestion in communication networks with arbitrary topology 
where all data sources use window flow control. Network 
elements such as queues and sources, are modeled as building 
blocks using the information conservation law which states 
that the information is either in transit, lost, or received. The 
proposed model implements the mechanisms ignored by the 
previously proposed models, notably at the queue and source 
levels. It is generic and independent of transport protocol 
specific congestion control algorithms. Previous models from 
the literature can be recovered from exact reduction or ap- 
proximation of this new model. The results obtained from the 
model match perfectly with the ones obtained from packet- 
level simulations. 

Future works will be devoted to the modeling of data loss, 
such as packet drops, and the time-out mechanism at the user 
level. 
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Figure 6: Predicted output flows. Plain: output flows, dashed: input flows 
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Figure 7: Model predicted and NS-2 simulated input and output number of 
packets: flow 1 (top) and flow 2 (bottom). 
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Figure 10: Scenario 1: Queue size 
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Figure 11: Scenario 2: Queue size 
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Figure 13: Scenario 3: queue 1 (top) and queue 2 (bottom) 
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Figure 14: Scenario 4: queue 1 (top) and queue 2 (bottom) 
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Figure 15: Scenario 5: queue 1 (top) and queue 2 (bottom) 
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Figure 16: Scenario 6: queue 1 (top) and queue 2 (bottom) 
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Figure 17: Scenario 7: Queue size (top) ACK buffer (bottom) 
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Figure 18: Scenario 8: Queue size (top) ACK buffer (bottom) 
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